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ABSTRACT 

We present a first estimate based on a cosmological gasdynamics simulation of galaxy cluster radio 
halo counts to be expected in forthcoming low-frequency radio surveys. Our estimate is based on 
a FLASH simulation of the ACDM model for which we have assigned radio power to clusters via a 
model that relates radio emissivity to cluster magnetic field strength, intracluster turbulence, and 
density. We vary several free parameters of this model and find that radio halo number counts vary 
by up to a factor of two for average magnetic fields ranging from 0.2 to 3.1 /iG. However, we predict 
significantly fewer low-frequency radio halos than expected from previous semi-analytic estimates, and 
we find that upcoming surveys will have difficulty distinguishing models because of large uncertainties 
and low number counts. Additionally, we find that expected number counts can be degenerate with 
both re-acceleration and hadronic secondary models of cosmic ray generation. We find that relations 
between radio power and mass and X-ray luminosity may be used to distinguish models, and by 
building mock radio sky maps we demonstrate that surveys such as LOFAR may have sufficient 
resolution and sensitivity to break this model degeneracy by imaging many individual clusters. 
Subject headings: cosmology:theory, structure formation, methods: N-body simulations, clusters of 
galaxies 



1. INTRODUCTION 

Although diffuse radio halos were discov ered in clus- 
ters of galaxies more than 50 years ago (Large et al. 
19591, complete statistical information about them has 



only been forthcoming within the past decade, owing to 
their rarity, steep spectra, and low surface brightn esses 



Radio survey s using the Very La r ge Array (VLA) ( Gio 



vannmi 



ner &: Sarazinli 



1999[ Cohen et al. 2007 Giovannini 2009 )71^6 



Westerbork Sy nthesis Radio Telescope (WSRT) ( Kemp 
m|2001], and the Giant Metrewave Radio 
MRTTi 



I'elescope (GMRT) ( [Brunetti et al.||2007[ ) have detected 
20 radio halos at redshifts up to z ~ 0.4, along with a 



variety of smal ler-scale radio features in clusters ( Kemp 
ner et aL|2004l ). Only about 1/3 of massive (> 10^^ Mq) 
clusters are known to host radio halos, and the halos 
themselves are not associated with any particular mem- 
ber galaxy, but rather dispersed througho ut the intra- 
cluster medium (ICM) ( [Feretti et aT]|20Ci4[ ). For clusters 
that do host halos, strong correla tions are see n between 
radio power and X-ray luminosity (|Liang et al.|200 0' Bac-, 
chi et al. 2003; Cassano et al. 2006; Brun etti et aL|2007p , 
halo mass (.Cassano et al. . 200 6 ) , and gas velocity disper- 
sion ( Cassano et al. |2008 ) . Also, observations indicate 



The proximate cause of diffuse radio halos is most likely 
synchrotron emission by high-energy electrons in galaxy 
cluster magnetic fields, but the means of generating and 
accelerating these electrons remains an open question, 
since these electrons have r elatively short (~ 0.1 Gyr) 
lifetimes. Dennison (11980') proposed that cosmic ray 
(CR) electrons are proHuced as secondary particles by 
collisions of > 1 GeV cosmic ray protons with ambient 
thermal ICM protons. The CR protons can be acceler- 
ated by shocks and quickly (< 1 Gyr) diffuse throughout 
the cluster; because of their larger mass, they have much 
longer synchrotron lifetimes than the electrons. This 
naturally ex plains the d i ffuse, cluster-wide properties of 
radio halos 
20 



1999 Pfrommer et al. 



mscuss a way in which 



2008| ). pCeshet 
he correlation 



&Loeb[ 

between radio and X-ray surface brightness can be ex- 
plained by hadronic secondary-type models, although 
this approach requires an extremely strong magnetic field 
(~ 5 ^G at radu ~ 1 Mpc for z = Q.2 - 0.4). Gamma- 
ray observations place limits on the abundance of high- 
energy protons, since in addition to producing charged 
pions that decay into the secondary electrons, the proton- 
proton collisions pro duce neutral pions, which decay into 



a str ong connection betwe en the presence o f a hal o and \^'^'^Y^t^^'!v photons ( |Blasi fc Colafrancescojl999| [Wolfe 



morphologi cal evidence for recent mergers (Buote 2001 
Brunetti et al. 2009J_jalthough some exceptions do~ex 
ist ( Russell et al.||201l] ). Indeed, recent simulations of 



et al. 20081. These observations indicate that hadronic 
cosmic rays may contribute at most 5-10 % of the to 



merging clusters suggest that the fraction of turbulent 
clusters is roug hly equal to the fra ction of clusters host- 
ing radio halos ( Vazza et al.||2011 1 . 
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tal pressure support in clusters (Ack ermann et al.|2010 ). 
However, the hadronic secondary model may have dif- 
ficulty explaining the diminishing of radio power with 
radius within the Coma cluster and the scaling of radi o 
power with X-ray luminosity (Donnert et al. 2010a|b ). 
Also, the ability of the high-energy protons to stream 
away from their sources can have significant implica- 



2 



tions for the resulting radio emission ( |EnBlin et al.|2011[ ). 
While gamma-ray observations have not strictly ruled 
out this model, recent estimates indicate a tension be- 
tween the magnetic field strengths r equired for this see 



nario and those observed in clusters ( Jeltema & Profumo 



Another promising explanation for the acceleration of 
the cosmic ray electrons is second-order Fermi accelera- 



tion by intracluster turbulence (Schlick cis_er et al. 1987 
Petros ian"2001| [Brunetti fc L azarian| |20iT| ). 'I'he cosmic 



rays themselve s must be injected into the ICM by radio 
galaxies ( JafFe] [T977 ) or accelerated by merger shocks. 
Since the lifetirne ot synchrotron-emitting GeV electrons 
in the intracluster ma gnetic field ((i?) ~ 1 — 10 ^G) 



is at most ~ 0.1 Gyr (Kuo et al. 2004 



20091, these electrons must be re-acce 



Brunetti et al 
erated by some 
process that operates in a more diffuse fashion. Hence a 
local acceleration mechanism is favored. The existence of 
some steep-spectrum low-frequency radio halos may sup- 
port the re-acceleration model, since the spectral index 
of the re- accel erated electrons dep ends on the turbulent 
Mach number ( Brunetti et al. 2008 ) , although future low- 



frequency surveys are necessary to fully test this picture. 
Assuming that ICM turbulence locally accelerates cos- 
mic rays to produce clusterwide radio halos, we expect 
that the radio emission should correlate spatially with 
the turbulent pressure (pw^). Indeed, this expectation 
is consisten t with observations using Chandra tempera- 
ture maps (Govoni et al. 2004). However, the fraction 



of CR electrons that are reaccelerated primaries rather 
than hadronic secondaries remains to be determined. 

Within the next decade and a half the number of 
known radio halos should increase dramatically owing 
to the development of sensitive low-frequency radio sur- 
veys based on technology originally developed for cel- 
lular communications. Examples of the operating and 
planned instruments include the LOw Frequency ARray 
(LOFAR), GMRT, the Karoo Array Telescope (KAT), 
and the Square Kilometer Array (SKA). LOFAR, for ex- 
ample, will be sensitive to radio frequencies between 20 
and 240 MHz and will be able to detect sources as faint as 
0.4 



2003 
ideal 



110 mJy at 15 - 



Rottgering et al ||2006[ ) 
or detecting and coun 



2 40 MH z ( Jarvis 
These c. 



2007 Rottgering 



laracteristics are 
ing radio halos such as the 



Coma radio halo ( ^640 mJy at 1.4 GHz with a spectral 
index ~ 0.4 — 0.75; Deiss et al.|1997j) as far away as a red- 
shift of 0.75. Moreover, if the CR electron population is 
dominated by reaccelerated primaries, then many more 
radio halos may be detected at low frequ encies owing to 
their steep spectra ( Cassano et al.pOO S*). 

Counts of cluster radio halos, in addition to probing 
the evolution of cluster merger activity, also potentially 
provide an additional means for using clusters to con- 
strain cosmological parameters. Unlike other methods 
for using clusters as cosmological probes that are based 
on their mass function or gas fraction, this measure is 
linked to their "instantaneous" formation activity rather 
than their time-integrated numbers. In principle its de- 
pendence on the cosmological volume element dV/ dz and 
the growth factor of linear density fluctuations Dj^{z) 
should also be different from and thus complementary 
to the more traditional measures. Additionally, if the 
cosmic-ray electrons responsible for the halos are acceler- 
ated by shocks and/or turbulence generated by mergers. 



it is reasonable to expect that recently merged clusters 
would display the most radio activity. Thus determining 
the abundances, spectral distributions, and other char- 
acteristics of radio halos as functions of redshift could 
provide information about the evolution of clusters and 
their merging activity over time. 

In this paper we present results from a numerical sim- 
ulation of cluster formation intended to study the form 
and evolution of the radio halo population as might 
be observed in a typical LOFAR survey. We apply a 
model of radio power that is generalized to include both 
hadronic and re-acceleration cosmic ray generation mech- 
anisms. We specifically choose this model to be as broad 
possible to allow immediate comparisons based on the 
same simulation. The means employed in the conversion 
of cluster density and velocity information into a simu- 
lated LOFAR radio sky are somewhat rudimentary given 
the uncertainties in the physics responsible for radio ha- 
los and the small scales on which it likely operates. How- 
ever, our results are the first based on combining a large- 
scale cosmological gasdynamics simulation with observed 
features of radio halos, and they show that future simu- 
lations with higher resolution and more realistic physics 
should enable straightforward comparisons with results 
from low-frequency radio observatories. While earlier an- 
alytical studies have involved more sophisticated models 
of CR generation, such as including spectral steepening 
effects, these have relied on the X-ray luminosity func- 
tion combined with the kno wn correlation between radi o 
power and X-ray luminosity (lEnfilin & Rottgering' "2002^, 
the Press-Schechter mass function (Cassan o et al.^2006) , 
or Monte-Carlo realizations ( jCas sano 2010|, whereas our 
work is based on di rect access to the intern al state of the 
simulated clusters. Donnert et al.| ( |2010a[ ) introduced a 
simulation of the Local Group including magnetic field 
injection and CR generation, but these results may be 
sensitive to the assumptions made about the magnetic 
field injection and are limited in volume. Our simu- 
lation, while involving simpler physics, covers a large 
(1 Gpc) volume, which will allow us to gather reli- 
able statistics and produce mock whole-sky radio maps, 
which are difficult to produce accurately with methods 
based on analytical mass functions. 

In Section[2]we provide details of the cosmological sim- 
ulation, while in Section[3]we explain the procedure used 
to associate a radio power with each cluster. We explore 
the range of valid models in Section |4] and use these re- 
sults to produce radio power relations in Section [5] Sec- 
tion [6] discusses the results in terms of radio halo counts 
as functions of flux and redshift, and we present various 
example radio sky maps in Section [7j Finally, we con- 
clude in Section [8] with a discussion of future directions. 

2. THE SIMULATION 

We simulated structure formation using the ACDM 
cosmological model within a periodic box spanning 
1024 h^^ comoving Mpc. We assumed a Hubble constant 
Hq = 100ft. km s"i Mpc"^ with h = 0.719, a present- 
day matter density parameter ri„j q = 0.262, baryonic 
density parameter J7{,,o = 0.0437, vacuum density pa- 
rameter r^A.o = 0.738, and spatially fl at geometry, as 
sugge sted by results from WMAP data (Komatsu et al. 
20111. Initial conditions for 1024^^ dark matter particles 
at a starting redshift Zi = 66 were generated using a 
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version of GRAFIC ( Bertschinger||200 l]) modif ied to ac- 
cept power spect ra generated by CMBFAST (Seljak & 
Zaldarriaga 1996). We normalized the power spectrum 
usmg (Tg = 0.74. We included adiabatic gasdynamics for 
the baryons using a perfect-gas equation of state with 
adiabatic index 7 = 5/3 and mean particle mass de- 
termined using interpolation from coUisional ionization 
equilibrium tabl es for primordial gas from [Sutherland 
& Dopita ( 1993 ) . Although no additional physics is in- 
cluded in this hrst calculation, we initialized the gas tem- 
perature at Zi to a constant value of 9100 K, correspond- 
ing to a preheating entropy of 250 keV cm^ at a redshift 
of 3. This level of preheating is adequate to reprod uce the 
observed X -ray luminosity-temperature relation (Bialek 
et al.||2001) for clusters of galaxies, although details of 



the scatter in this relation and its correlation with other 
cluster properties such as the presence of cold cores are 
not constrained to match observations. 

We ran our simulation using the FLASH code version 
3.3 dFryxell et aL| |2000'; Dubcy ct al. '2008^ using a new 
direct multigrid Poisson solver (Ricker 2008| ) with 1024^^ 
dark matter particles and a unif orm 1024'^ base mesh. 



The piecewise- parabolic method (Colella &: Woodward 
19841 was used to solve the Euler equations ot gasdy 



namics. To achieve the resolution necessary to estimate 
the level of turbulence within clusters, we used adaptive 
mesh refinement within 100 preselected regions. Each 
region was 50 Mpc on a side centered on a halo 
identified using a lower-resolution precursor run. The 
halos were selected to uniformly sample the range of re- 
solvable halos using mass function weighting to ensure 
a representative sample in the full simulation. Within 
the preselected regions we used a dark matter particle 
refinement criterion, allowing no more than 100 dark 
matter particles within a zone. We refined to a maxi- 
mum resolution of 32 kpc. Estimates of the integral 
scale of turbulence in clusters suggests a power law spec- 
trum from spatial scales of 0.8 to 8 kpc, with no visible 
turnover (Kuchar & Enfilin 2011), suggesting that the 
integral scale is larger than 8 kpc and that our resolu- 
tion should allow us to determine the level of turbulent 
dissipation. Since the refined regions were larger than 
the halos on which they were centered, we captured a 
total of 131 high-resolution clusters. We ran the code 
on the Cray XT5 machine at Oak Ridge National Lab- 
oratory, where the simulation required approximately 
450,000 CPU-hours on 16,000 processors. Output files 
containing both particle and gas information were writ- 
ten beginning at z = 2.0 at every Az = 0.25 for the 
purposes of mock sky generation. 

2.1. Halo finding and virial mass 

We created halo catalogs from the simulation outputs 
using the friends-of-friends (EOF) technique with a link- 
ing length parameter b = 0.2 and considered only halos 
with at least 3, 000 particles (i.e. an FOE dark-matter 
mass of 2 X 10^"* h^^ M0). Our base-grid spatial resolu- 
tion is sufficient to ensure accurate counts of halos with 



consider here (Heitmann et al. 2005 



Figure [T] shows our mass function for all halos in t 



Lukic et al. 


2007 


all halos in tJ 


le Sim 



ulation volu me compared against the best fit of [Warren] 
et al.| ( |2006[ ). We find fewer high -mass objects relative 
to the Warren fit, but this is not unexpected (see, for 




15.5 



Log [ M| 



Fig. 1. — Mass function of all halos in the si mulation volume 
compared against best fit of [Warren et al.] | |2006[l . Errors bars are 
given at 2a and the vertical arrow denotes our FUP 
ness limit. 



JF halo complete- 



example, Knebe et al. 2011). Also, we tend to over- 
produce low mass objects, even below our resolvability 
limit. However, note that even though we produce too 
many low-mass objects relative to the Warren fit below 
our resolution limit, our mass function still turns away 
from the expected slope, and thus we cannot fully trust 
the number counts below this threshold. There were 
^ 4000 resolvable objects at z = 0.0. 

To make comparison s with the observational analysis 
of Cassano et al. (2006) (hereafter CBS06), we compute a 
spherical overdensity radius Ry for each of our halos. For 
the high-resolution sample (i.e., the adaptively- refined 
halos within the 100 predefined regions), we compute 
overdensities including both just dark matter and with 
dark matter plus gas. For the remaining fixed-resolution 
halos outside the predefined regions, we only include dark 
matter in the overdensity calculation since the gas data in 
these halos were poorly resolved. In Section [6] we discuss 
our procedure for assigning radio power to these lower- 
resolution halos. We use the same definition of ove rden- 
sity as in CBS06, namely Kitayama fc Suto| ( [1996 ): 

(1) 



Ae = 187r2 (1 -I- 0.4093w(z)0-9"52) ^ 
where cj(z) = 17^ (z)^^ — 1. Here, 

The virial mass. My, follows as My = (4/3)7rAcPm (z)i?^, 
where Pm{z) is the mean mass density: 

p,„(z) = 2.87 X 10"0^,o(l + zf Mq Mpc'^. (3) 

The most massive cluster in our simulation has a mass 
My = 1.2 X 10^^ Mq. 

3. SIMULATING RADIO EMISSION 

We identify gas zones within Ry for each halo and 
associate them with the halo in which they are found. 
We create two-dimensional maps of projected density 
and projected turbulent pressure, niivf, where is the 
mass in the cell i and the average velocity is defined 
as the difference between the measured velocity in the 
cell and the bulk velocity averaged over 300 kpc regions, 
Vi = Vi — V300 kpc- This is a strategy for removing bulk 
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motions similar to 
proach to re move 



Nazza, et al. (2011 ) (an alternative ap- 
3ulk velocities from this calculation is 



discussed in Paul et al. 2011). We use these projections 



in two ways: to create simulated surface brightness maps 
and to construct total radio luminosities by integrating 
these quantities across the entire projected cluster sur- 
face out to the virial radius. The integrated projected 
density is of course My and we will designate the inte- 
grated turbulent pressure as Ty = J2i "t-i^^I- 

Therefore we can create simulated radio surface bright- 
ness maps for our clusters by normalizing maps of pro- 
jected turbulent pressure and projected mass using as- 
sumed radio luminosities and rest-frame spectra. Be- 
cause we may not fully resolve intracluster turbulence, 
the total amount of turbulent pressure in our clusters 



seen in high-resolution 


simulations dRicker & Sarazin 


2001 [Ritchie & Thomas 


2002| Mitchell etai.|2009|). How- 


ever, because the norma 


lization of the radio power is sup- 



and Bqmb = 3.2(1 



(5) 



plied independently (see below) , all we require of the tur- 
bulent pressure maps is that they be sufficiently diffuse 
and representative in spatial extent of clusters containing 
radio halos. The detailed structure of the maps should 
not be regarded as realistic. Because of the beam smear- 
ing described below, this fact does not significantly affect 
our analysis. Also, while we cannot depend on these sim- 
ulated clusters to provide correct high-resolution X-ray 
and radio surface brightness maps, we can still use them 
to identify broad features, such large-scale shocks and 
the relative radial dependence of turbulence. 

Using these integrated quantities, we construct a rest 
frame 1.4 GHz radio power via 

PiA GHz = C.B.M^T'y, (4) 

where Cs is a scaling constant. My is the virial mass, r^, 
is the virial turbulent pressure, and Bg is the magnetic 
field parameter: 

' {B{MyY 4 
where B{My) = (B) {My / (M))'' 

2)^ fiG is the equivalent magnetic field strength of the 
cosmic microwave background. This formulation sepa- 
rates physical processes that generate cosmic rays {My 
and r„) from those that contribute to radio emission 
{B in the numerator) and CR losses due to emission 
{B^ in the denominator) and inverse Compton scatter- 
ing (iJcMs)- The losses enter into this equation because 
they limit the maximum CR energy. In this formal- 
ism. My measures the total cluster mass and thus should 
scale with the dependence of cosmic ray generation on 
hadronic secondary processes, whereas Ty measures the 
total cluster turbulence and thus should provide a mea- 
sure of the re-acceleration of cosmic ray electrons by that 
turbulence. We will set (M) = 1.5 x 10^^^ Mq. There are 
thus five independent parameters: the average magnetic 
field (B), the scaling of magnetic field with cluster mass, 
6, the scaling of radio power with virial mass, a, the scal- 
ing of radio power with turbulent pressure, c, and an 
overall scaling parameter Cg ■ A summary of our model 
parameters is given in Table [ij 



jCassano & Brunetti 


( 


2005 


1. Note that the 


ironic secondary models of 


Dolag & Enfilin 



TABLE 1 

Parameters of the radio luminosity model. 



Parameter 



Description 



Cs Overall scaling 

(B) Average magnetic field 

b Scaling of magnetic field with cluster mass 

a Scaling of radio power with cluster mass 

c Scaling of radio power with turbulent pressure 



E 0.5 : 




14.4 14.6 14.e 
Log [M„/ Mel 

Fig. 2. — Total turbulent pressure, F^, versus virial mass, My 
at z = for the high-resolution cluster sample. Shown is a best-fit 
line in log space. The slope of the line is ^ 1.7. 



( 2000 1 identified the functional form of the magnetic field 



i.e. the denomina- 



as~B^ B{MyYI{B{My) 

tor is not squared. Our model easily accommodates this 
scenario: when Scmb dominates, this will appear as a 
constant factor folded into Cs, and when B{My) domi- 
nates this will simply adjust the mass scaling factor a. 
We are fixing the form of the magnetic field dependence 
since the radio synchrotron power will always depend on 
magnetic field pressure {B^) independently of t he cosmic 
ray generation and acceleration mechanisms (see Cassano 
& Brunetti (2005) for a discussion). We stress that al- 
though this model is relatively simple a nd ignores the 
stream ing and diffusion of cosmic rays (Enfilin et al. 
2011), it allows us to explore a range of plausible ac- 
celeration mechanisms and examine relative changes to 
luminosity functions, scaling relations, and other radio 
properties. 

This model allows us to explore both cosmic ray gener- 
ation mechanisms or a mixture of both. For example, the 
hadronic secondaries model should predict radio power 
which scales with cluster mass, so c = in this case. A 
re-acceleration model is proportional to turbulent pres- 
sure, so a would be 0. Note that the model of CSB06 
is based on re-acceleration, but only scales with cluster 
mass. This is because Vy roughly scales with My with 
a logarithmic slope of 1.7, as shown in Figure [2] There 
is some scatter in this relation due to the merger history 
of a particular cluster: recent mergers produce stronger 
levels of turbulence, which tend to scatter the cluster 
higher in this relation. Note that the model of CSB06 
corresponds here to a = 4/3 and c = 0. 

A degeneracy exists for calculations of total radio lu- 
minosity between models that scale with turbulent pres- 
sure and those that scale with mass, since we may freely 
exchange c for 1.7a and vice- versa. However, a more de- 
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tailed examination of cluster atmospheres reveals striking 
differences. Even with the relatively low resolution of our 
simulation, and the resulting inability to fully reproduce 
correct structures in the cluster atmospheres, we can 
identify gross differences in the projected maps. Figure 
[3] shows projections of mass and turbulent pressure for 
two clusters. The mass projections of both clusters are 
roughly spherical, as expected. However, the turbulent 
pressure maps show more varying morphology. While the 
halos are roughly equal in mass (~ 8 x 10^'' h"^ M©), 
one shows much greater turbulent structure, indicating 
recent merger activity, which may explain the scatter 
in Figure [2] Thus, even though different mechanisms 
of CR generation may produce similar cluster counts (as 
we will see below) , high-resolution radio and X-ray imag- 
ing of clusters may help to determine which mechanism 
dominates. 

Since we do not include in our simulations any de- 
tailed cosmic ray generation mechanisms, and we because 
we want to keep our model as general as possible, we 
must fix the scaling parameter Cg by using observations. 
This scaling will then combine any extra constants and 
parameters not included in our analysis. For a given 
set of model parameters, we set Cs by assigning a ra- 
dio luminosity to the most massive cluster in our sim- 
ulation. This cluster has a mass '-^ 3 x 10^^ M©, 
which fits within the observed radio halo mass range of 
2x1015-6x101^ H-'^Mq. We do this with the Pi. 4 -M^ 
relation found in CBS06, which is based on combining the 
observed correlation of radio halo power and X-ray lumi- 
nosity with the correlation between X-ray luminosity and 



log 

(2.9 ±0.4) log 



M„ 



1015 /i-i Mq 



(0.814 ±0.147) 



(6) 



We then apply this same constant scaling to all remain- 
ing high-resolution halos in the sample. While the scal- 
ing may contain some additional dependence on mass or 
turbulent pressure not accounted for in our parameter- 
ization, this can easily be accommodated in our study 
by adding to (or subtracting from) the parameters a and 
c. An example of a particular model compared against 
the observed relation is shown in Figure |4] Our best-fit 
relation in this plot and throughout this paper uses only 
the clusters within the observed mass range. 

Since no radio halos have been observed beyond a red- 
shift of ~ 0.4, and available statistics do not strongly 
constrain evolution in this relation, we will fix the scal- 
ing at z = and apply the same scaling to higher-redshift 
clusters. We will also assume power-law energy spectra 
with a spectr al index of 1.2, cons istent with low-redshift 
observations ( jFeretti et al.||2004| . Finally, we do not in- 
clude in our niodel the relationship between synchrotron 
break frequency and the presence of a radio halo, which 
can be used to calibrate models to the observed fraction 
of clusters hosting radio halos (CBS06). We will discuss 
the potential impacts of this assumption in the conclu- 
sion. 

4. EXPLORATION OF VALID MODELS 

To constrain our model choices we make selections for 
the model parameters, assign radio powers to the clus- 
ters using the procedure described above, find the best fit 



line to our derived Pi. 4 — data above a mass thresh- 
old of 10^5 Mq, and compare the best-fit slope and 
normalization to the observed values. We only accept 
model choices that produce fits that lie within la of the 
observed relation. This is a strategy similar to the one 
employed by CBS06, but that we are applying a slightly 
more stringent test by enforcing the known relation to 
somewhat lower radio powers than they consider. We do 
this so that we can capture enough halos 10) to gen- 
erate sufficient statistics for our best-fit lines. Obviously, 
we could just select two models that span the valid range 
and analyze their difference, but we wish to explore the 
relationships among the various model parameters and 
the separate consequences of varying each one. 

Figure [5] shows colormap plots of allowable models. We 
vary (P) from 0.2 to 6.0 ^G, b from 0.5 to 1.5, a from 
0.0 to 5.0, and finally c from 0.0 to 3.0. We could ex- 
plore even larger values of a and c, but as we will discuss 
below la uncertainties in the measured Pi.4ghz — 
relation place upper limits on the scaling of a and c at 
these chosen maximum values. We also assume a positive 
correlation between radio power and My and r^,. While 
we allow the mass and turbulent pressure scaling param- 
eters to vary all the way to 0, we constrain the scalings 
associated with magnetic fields. We constrain the aver- 
age cluster magnetic field strength from 0.2 /iG, which 
is set by observed upper limits on hard X-ray emission 
(CBS06), to 6.0 /xG, which is a reasonab le upper limit 



(e.g. 



Govoni et al. 



from rotation measure obs ervations 
2004| [Bonafede et al.| | 201 1 [). The rest r iction s on b come 



from the simulations of Dblag et al. (2002), which fol- 
lowed the adiabatic compression of seed magnetic fields 
as clusters formed. They found a scaling P oc M^-^^. We 
allow some uncertainty in this value, but do not allow 
a compete lack of scaling of magnetic field with cluster 
mass. For simplicity, we have combined the mass and 
turbulent pressure values as a ± c, so a ± c is varied from 
0.0 to 9.0. The contours for each individual parameter 
show structures similar to those for this combined pa- 
rameter. In these plots we are showing the maximum 
allowed value for a given point on each contour plot. All 
values less than the plotted value are also allowed. For 
example, for b = 1.0 and (P) = 3.0 /iG the allowable 
values for a + c are from 0.0 to 1.5. 

We find that very strong magnetic fields are only al- 
lowed if the scalings with virial mass and turbulent pres- 
sure are very steep. In these cases strong radio power in 
low mass objects due to high (P) is offset by significantly 
lower radio power associated with or F^ . If the scal- 
ing of magnetic field strength with cluster mass is above 
unity, then it is difficult to fit strong magnetic fields at 
high mass within the observed relations. We find several 
regions forbidden in our models: strong magnetic fields 
coupled with low a + c, and very low or very high a + c 
and b values. 

We see interesting structures in the contours: steps 
and wiggles in the a + c plots, and striations in the oth- 
ers. These are due to the scatter that develops in the 
Pi. 4 — My relations and the resulting variations of the 
best fit lines. Because of this variation, we do not see 
monotonically increasing (or decreasing) behavior in the 
contour plots, especially at extreme values. Surprisingly, 
we find that a = c = 0.0 is allowed, but only at low 
(P) and high b. This is because of the implicit mass de- 
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Fig. 3. — Projected density (left-hand plots) and projected turbulent pressure (right-hand plots) for a two clusters (top and bottom 
rows). Projections are taken along the z-direction within Rv for each cluster. The units for projected density are 10** g cm~^ and for 
projected turbulent pressure are images are normalized to a uniform grid 50 cells on a side. The corresponding comoving 

length scale is indicated in the figure. 

pendence in the calculation of the cluster magnetic field 
strength, B{My). Also, the model used in CBS06 is for- 
bidden in our analysis, since we are enforcing the known 
relation to a wider range of cluster masses. However, we 
can allow their model choice when restricting ourselves 
to the mass ranges they consider. 

We use these contours to guide our selection of models 
for further study. We wish to adequately sample the 
space of allowable models and explore the limits allowed 
by observational constraints. We also wish to explore the 
effects of holding one parameter constant and varying the 
others to their extreme allowed values. To aid analysis, 
we collect our choices into six model groups, enumerated 
in Table [21 In this table we list the values chosen for 
a particular parameter set and a unique designation for 
that set used in further plots. 

In Model Group 1 we set a to 0.0, fix c = 0.7, and 



vary the magnetic field parameters as widely as possi- 
ble from a minimum of {B) = 0.2 to 3.1 /iG. We also 
vary the scaling parameter associated with the magnetic 
field, h. Each {B) is coupled with a unique 6, except for 
{B) = 1.5 /iG, where we examine h = 0.8 and b = 1.1, 
which are the minimum and maximum allowed values 
for this particular configuration. We explore the oppo- 
site behavior in Model Set 2 by fixing the magnetic field 
parameters to {B) = 2.0 and simple linear scaling 
6 = 1.0 while having no explicit dependence and 
studying the maximum and minimum allowed values for 
turbulent pressure scaling, c. We chose this value of the 
magnetic field so that we could get the maximum differ- 
ence in c. We keep the same magnetic field configuration 
for Model Set 3, but here we exchange a and c using the 
measured relation (Figure [2]) . In this Model Set we fix 
the quantity a -f 2c. This allows us to hold the magnetic 
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Fig. 4. — Radio halo luminosity versus virial mass for one exam- 
ple model (points) with best fit above 10^^ Mq (solid line) 
compared against the observed best fit found in CBS06 (thick 
dashed line) and Icr uncertainties (thin dashed lines). Points are 
the high-resolution cluster sample. The example uses parameters 
(B) = 2.0 AtG, b = 1.0, a = 0.0, and c = 0.7. 



TABLE 2 
Model groups and parameter sets. 



Designation {B){fiG) b a c 



Model Group 1) Fixed c. Varying Magnetic Field 



lA 


0.2 


0.650 


0.000 


0.700 


IB 


0.5 


1.000 


0.000 


0.700 


IC 


1.5 


0.800 


0.000 


0.700 


ID 


1.5 


1.100 


0.000 


0.700 


IE 


3.1 


1.470 


0.000 


0.700 



Model Group 2) Fixed Magnetic Field 
2A 2l) 1.000 0.000 0.510 

2B 2.0 1.000 0.000 0.930 

Model Group 3) Exchanging a and c 
3A ■U) 1.000 1.300 0.000 

3B 2.0 1.000 0.650 0.325 

30 2.0 1.000 0.000 0.650 

Model Group 4) Extreme Allowed Magnetic Fields, Fixed c 
4A 02 1.000 0.000 0.700 

4B 3.1 1.000 0.000 0.700 

Model Group 5) Extreme Allowed Magnetic Fields 
5A 02 0.500 0.000 0.875 

5B 6.0 1.500 2.000 0.000 

Model Group 6) Extreme Allowed Scalings 
6A 02 1.260 0.000 0.010 
6B 6.0 0.600 0.000 1.550 

field fixed while going from a hadronic-like CR model 
(c = 0) to a re- acceleration model (a = 0). We designed 
this Model Set to verify that our results are robust to 
even exchanges of a and c using the measured relation, 
which they should be. In Model Set 4 we fix 6 = 1.0, 
a = 0.0, and c = 0.7 and examine the extreme allowed 
average magnetic field. We chose these values of b, a, 
and c such that we could get the maximum change in 
(B). We repeat this test in Model Set 5, but now allow 
6, a, and c to vary to accommodate the extreme val- 




0.6 0.8 1.0 1.2 1,4 

b 

(a) a + c as a function of (B) and b. 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 

a + c 

(b) (£?) as a function of b and a + c. 




1 2 3 4 5 6 

<B> 

(c) 6 as a function of a -f c and (B). 



Fig. 5. — Contours of allowed radio halo model parameters. Al- 
lowed models are determined by fitting a line to our Pi. 4 — M„ data 
above M^^ Mq and ensuring that the slope and normalization 
are within la of the known relation. {B) is given in units of /iG. 
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ues studied of (B). Finally in Model Set 6 we pick two 
model parameter sets that represent extremes of all four 
parameters. 

5. RADIO POWER RELATIONS 

We begin our analysis by using our sample of 131 
high-resolution clusters to examine the relationship be- 
tween radio luminosity and virial mass and X-ray lumi- 
nosity. To construct X-ray luminosities, we use the mekal 
plasma emissivity model supplied with the XSPEC pack- 
age (Arnaud 1996). We then build a composite X-ray 
spectrum for each cluster and use that spectrum to gen- 
erate the 0.1 — 2.4 keV rest-frame luminosity within the 
spherical radius for each cluster. Since our simulation 
does not include cooling and central AGN feedback, our 
cluster temperatures, and hence X- ray luminosities, are 
uniformly higher than observed (see Stanek et al. (20101 
for a discussion of such effects). However, the slopes 
of our relations are still within observed limits, and we 
can still study the relative differences among models and 
their evolution with redshift. 

In Figure |6] we show the best-fit slope and normaliza- 
tion for each model, grouped by model group, for the 
Pi A — My relation generalized from (Eq. rol): 



log 



Pi. 



3.16 X 1024/i-p^W Hz" 



Af log 



1015 hjo' M, 



'-0 



(7) 

where Af and bf (note that we have added the sufhx / 
for "fit" to distinguish these from the parameters used 
in our radio power model) are the slope and normaliza- 
tion, respectively. We fix the scaling parameter Cg at 
each redshift for each model. In essence this assumes 
that the Pi.4ghz — My relation holds even at high red- 
shift. For each model we show three points: one each for 
z = 0.0, 0.25, and 0.5. Above redshift 0.5 we do not have 
enough halos above the minimum mass threshold to gen- 
erate meaningful statistics. Note that the observational 
uncertainties essentially fill the entire plotting space. 

For Model Group 1, where we vary only the magnetic 
field parameters, we see a progressive steepening of the 
slope with higher b values, as expected. Note especially 
the differences between models IC and ID, which have 
identical values of (B). The values of b and (B) also 
jointly affect the normalization of the -Pi.4ghz ~ My re- 
lation, with smaller values b generally leading to lower 
normalizations. The redshift evolution of the models in 
this Model Group shows diverse behavior. For models 
with b < 1.0 the normalization tends decreases with in- 
creasing redshift, but not very significantly. This makes 
sense as the clusters are in general uniformly smaller at 
higher redshift. However, the slope increases a.t z — 0.25, 
which perhaps suggests greater variance in the turbulent 
properties of the clusters. The uncertainties in the values 
for redshifts 0.25 and 0.5 make them difficult to compare 
against each other. With b > 1.0, the trend reverses and 
the normalization tends to increase with high redshift. 

In Model Group 2, where we keep the magnetic field 
fixed and vary the scaling with turbulent pressure, we 
see that, as expected, larger values of c lead to steeper 
slopes and lower normalizations in the best fit relation. 
Since our scatter is related to the turbulent pressure, 
models with higher values of c will have correspondingly 
larger uncertainties. We see similar redshift dependence 



for model set 2 A as in Model Group 1, but the model set 
2B displays reverse behavior (increasing normalization 
with redshift), although the uncertainties are so large as 
to make firm statements difficult. The steep dependence 
on turbulent pressure overwhelms the general mass de- 
pendence, so that even though the clusters are in general 
smaller at higher redshift (leading to a lower normaliza- 
tion), the turbulence in the most massive cluster (which 
is used to fix the normalization) increases, negating the 
mass effect. 

For Model Group 3 we see that exchanging a for c in 
the radio power model does lead to small differences in 
the slope and normalization. When we set c = (so that 
there is no dependence on turbulent pressure) we essen- 
tially eliminate the uncertainties. The small differences 
are due to the fact that our scaling relation between F^ 
and My is only a best-fit approximation, and that scat- 
ter in that relation can affect the resulting Piaghz ~ My 
relation. All the points, however, are within 2a of each 
other. All these models show identical redshift depen- 
dence. 

We see the drastic effects of changing the assumed av- 
erage magnetic field in Model Group 4. Higher magnetic 
fields lead to higher normalizations and flatter slopes. 
However, the redshift evolutions exhibit opposite trends, 
such that at z = 0.5 the effects of the magnetic field are 
indistinguishable from each other. For weak magnetic 
fields, the cosmic rays are dominated by their interac- 
tions with the CMB, and the equivalent pressure of the 
CMB increases with redshift, lowering the synchrotron 
power at higher z. For strong fields, the dependence on 
cluster mass is more explicit, and at higher redshifts the 
clusters are, in general, smaller and the mass distribution 
has a steeper slope. 

In Model Group 5 we see that the dependence on the 
scalings overwhelms the dependence on the average mag- 
netic field. Even though model set 5A has the same (B) 
as set 4A, the dependence on b, a, and c forces a much 
fiattcr slope. Thus we may conclude that the effects of 
average magnetic field are degenerate with the scaling 
parameters, although the parameters taken individually 
can lead to significant differences. Similar behaviors are 
displayed by Model Group 6, although the error bars are 
so small because of the weak dependence on the turbu- 
lent pressure. 

In Figure [T] we repeat the above analysis for the Pi 4 — 
Lr relation: 



log 



Pi.. 



3.16 X 1024 /i-iW Hz"^ 



Af log 



ijQ ergs s 
(8) 

Wc find similar behaviors as in the Piaghz ~ My rela- 
tion plots above, except that our error bars are generally 
larger due to scatter from our estimates of L^. Since X- 
ray luminosities are generally easier than virial masses to 
compute from observations, observations of many more 
radio halos may reduce the statistical uncertainties to 
such a level as to potentially distinguish the allowed scal- 
ings and dependencies. The redshift evolution of the 
Pi A — Lx relation especially may provide a way of de- 
termining the dominant components of radio power and 
the average magnetic strength of clusters, as we have 
discussed above for the Pi. 4 ~ My relation.. 
In general, if the observational uncertainties in A/ and 
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Fig. 6. — Best fits for the Pi. 4 — relation for each radio halo parameter set. Each best fit to Eq.[7| generates a slope and normalization, 
which we represent as a point with Ifj error bars. Solid lines are the best fit at 5; = 0.0, thick dashed lines are z = 0.25, and thin dashed 
lines are z = 0.5. We have identified each model with its designation from Table[2]and the portions of the model that change in the given 
model group. (B) is given in units of ^G. 
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Fig. 7. — Best fits for the Pi. 4 — Lx relation for each radio halo parameter set. Methods and colors are identical to Figure[6l We have 
identified each model with its designation from Table [2] and the portions of the model that change in the given model group. (B) is given 
in units of fiG. 



bf are reduced by approximately a factor of two, many 
degeneracies in the model parameters will be eliminated. 
We note that we are basing this analysis on our sample 
of only 131 clusters. While this is significantly more than 
the current known number of radio halos, it is still far 
fewer than we expect to see with instruments such as 
LOFAR, as we will see below. While more objects could 
reduce the uncertainty, the precise amount of error also 
depends on the intrinsic scatter in the observed relations, 
which can be affected by biases. We also assume a power- 
law relationship holds between radio power and cluster 
mass even to low-mass clusters, which may not be the 
case once significant numbers of low-luminosity halos are 



detected. 



6. LUMINOSITY FUNCTIONS AND RADIO HALO COUNTS 

We now turn to a discussion of these models in terms 
of total counts of all objects in the simulation. To do 
this, we must assign a rest-frame radio luminosity to each 
cluster in the simulation box, even if it is not at high res- 
olution. We accomplish this by combining the derived 
Ty — My relations described above with a relationship 
between M„, which includes gas, and Myj^yi, which only 
includes dark matter. Our smallest clusters are only a 
few zones across and thus do not contain enough gas 
zones to accurately capture the contribution of the gas 
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Fig. 8. — Data for (points) and best fit to (solid line) the My - 
M„ oM relation. Points are the high-resolution cluster sample. 

to My. However, we defined our FOF halo completeness 
limit so that we can always get reliable evaluations of 
My^UM (i-e., Ry > Ax, where Ax is the resolution of our 
pre-refinement uniform grid). This gives us an interpo- 
lated value of for each fixed-resolution cluster, which 
we then feed into Eq. |4]to generate a radio power for that 
cluster. 

We find a very tight correlation between My and 
My^DM for our high-resolution sample, as shown in Figure 
m We fit a line to these data and found the correlation 
to be 



log 
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IQis h 
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0.99 log 
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1015 h-^ 



Mr. 



+ 0.08. 



(9) 

This relationship implies a unifor m gas fraction consis 



tent with other simulations (e.g., Stanek et al. 2010). 
We use this fit to extract an My for each fixed-resolution 
cluster (i.e., those outside the refinement regions) which 
is then used to compute its equivalent radio power using 
the — My relation. While equivalent to directly inter- 
polating from a F^ — M^^dm relation, we found that this 
procedure produces less scatter and hence more reliable 
interpolations. We do not add additional scatter to the 
interpolation procedure. For the analysis below we will 
include Poissonian uncertainty where appropriate, and 
this dwarfs any uncertainty introduced by scatter. By 
binning our data for luminosity functions, the main ef- 
fect of scatter is to simply move clusters around within a 
given luminosity bin, and any clusters that are scattered 
into a luminosity bin are roughly offset by clusters scat- 
tered out of the same bin, especially at low luminosities 
where the function is relatively flat. 

Not every cluster hosts a radio halo, and for simplic- 
ity we will only assign radio halos to 5% of our clusters. 
This is roughly in line with observations, which indicate 
that 1/3 of clusters above 2 x lO^^ M(7) and only 
2 — 5% of smaller clusters host radio halos ( Cassano et al 



2008[ ). We do this rather than employing a radio halo 
probability that is a step function in mass because we 
do not have enough high-mass (> lO^^ Mq) halos to 
calibrate our number counts based on the intrinsic scat- 
ter in our derived Piaghz — My relation. Note that in 
the context of hadronic secondary models, this enforced 
fraction implies that not every cluster is sufficiently mag- 
netized to generate radio emission. More sophisticated 



techniques to calibrate number counts exist, such as us- 
ing the synchrotron break frequency, Vb , which is used in 
the identification of radio halos (see CBS06 for a discus- 
sion). However, again we do not have enough high-mass 
objects to use this approach. With enough high-mass ha- 
los, the precise calibration would depend on our choice 
of parameters (B), 6, a, and c. Thus we may be over- 
counting the number of radio halos at 1.4 GHz. However, 
we can still gather useful results as to the relative effects 
of varying radio halo models. Also, our analysis will in- 
clude results at 150 MHz, where Vf, is much lower and 
hence our results are more valid. To calculate the radio 
luminosity at 150 MHz, we assume a simple power law 
with spectral index 1.2. Our smallest resolvable cluster 
has a 1.4 GHz radio luminosity of 2 x 10^^ W Hz^^. 

Figure [9] shows our calculated radio halo luminosity 
functions at redshift z = 0.0 at 1.4 GHz and 150 MHz. 
At high masses our luminosity function is well be low the 
estimated values of Enfilin & Rottgering ( 2002 ) . How- 
ever, this is not unexpected due to our limited simula- 
tion volume and cosmic variance, and the fact that they 
assume a Press-Schechter mass function with a fixed 1 /3 
fraction of clusters hosting radio halos. For Model Group 
1, in which we fix the dependence on cluster mass and 
turbulent pressure but vary the average magnetic field 
strength and scaling of magnetic field with cluster mass, 
we see a bifurcation in the luminosity function at low lu- 
minosities: models with 6 < 1.0 produce up to a factor of 
two more low-luminosity radio halos than those models 
with b > 1.0. The distinction is much more significant 
at 150 MHz, where more objects allow for smaller uncer- 
tainties (assuming perfect detector sensitivity — a point 
which we will address later in the discussion of number 
counts). However, despite difference in the Pi.4ghz — My 
relations, models lA and IC are largely indistinguishable 
from each other, as are models IB, ID, and IE. Number 
counts at low luminosity are reduced significantly when 
the magnetic field scaling b is greater than unity. 

We see in Model Group 2, in which we keep the mag- 
netic field values fixed, that steeper scalings with tur- 
bulent pressure produce many fewer radio halos at both 
1.4 GHz and 150 MHz. The models at 1.4 GHz are only 
statistically distinguishable at the lowest luminosities, 
but the models at 150 MHz are easily separable through- 
out almost the entire range of radio luminosities. Similar 
conclusions can be made regarding Model Group 4, which 
brackets the extreme allowed magnetic fields with fixed 
scalings, where strong differences in the assumed average 
magnetic field strength lead to somewhat distinguishable 
differences. These behaviors persist for Model Groups 5 
and 6, in which all parameters are allowed to vary, with 
the general rule being that it is difficult to separate these 
models when only relying on 1.4 GHz halo counts. The 
only models that remain inseparable are those in Model 
Group 3, in which our exchanges of a and c lead to degen- 
erate radio halo counts at both frequencies, as expected. 

To accumulate total and binned counts we set a van- 
tage point in the center of our computational domain and 
find halos whose locations lie on the light cone emanating 
from this position. To do this we use saved checkpoint 
files at z = 0.0, 0.25, 0.5, 0.75, and 1.0. Moving out- 
wards in small (Az = 0.05) redshift slices covering the 
same range, we locate the nearest output (in redshift) 
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Fig. 9. — Radio halo luminosity functions at z = 0.0. Solid lines are luminosity functions at 1.4 GHz and dashed lines are at 150 MHz. 
Error bars indicate la Poisson uncertainties. We have identified each model with its designation from Table [2] and the portions of the 
model that change in the given model group. (B) is given in units of ^G. 



to each slice and use the center-of-mass peculiar veloci- 
ties of the clusters found in that output to estimate new 
positions in the slice under consideration. We then com- 
pute the flux as P/(47rc?|), where rfi, is the luminosity 
distance of the cluster. We do this at both 1.4 GHz and 
150 MHz assuming a spectral index of 1.2. 
We begin with Figure TOl where we show the total 



counts of radio halos at 1.4 GHz and 150 MHz in the 
observable universe as a function of flux limit in mJy. 
To generate error bars we propagate the la uncertain- 
ties in the derived M„ - M„^dm and r„ — My fits to gen- 
erate a minimum and maximum radio power for each 
cluster. Our number counts are bound by the resolv- 
ability limit of our simulation. While the model trends 
continue from the above analysis, we find that at high 
flux limits (> 100 mJy) and high frequencies, we have 
too few radio halos to strongly distinguish several mod- 
els, even those with large discrepancies in either assumed 
average magnetic field or scalings with virial mass or to- 
tal turbulent pressure. This is due to the suppression of 
radio halos at high redshift, meaning that the integrated 
counts depend most strongly on high-luminosity objects, 
where the counts are nearly the same. At 150 MHz and 
an assumed LOFAR sensitivity limit of 30 mJy, we find 
that although some models, such as model sets 2A and 
2B, produce an almost factor of two difference in the 
total counts, the large uncertainties preclude any clean 
distinction. 

In Figure [TT] we show the total counts of radio halos 
within redshift z < 0.2 at 1.4 GHz and 150 MHz. This 
redshift range fits largely within our computational vol- 
ume without the need for periodic replication of the do- 
main and is more easily accessible to observers. Although 
we find little degradation in the total number counts in 
the LOFAR-accessible regime (> 30 mJy), the models 
remain indistinguishable. 



7. SIMULATED RADIO SKY MAPS 



While we could in principle produce mock sky maps 
within any frequency range, we choose LOFAR-like pa- 
rameters since low-frequency instruments are able to sur- 
vey large portions of the sky and hence collect many halo 
images for use in statistical comparison. We generate 
raw mock sky maps in the 20 - 240 MHz LOFAR band- 
pass by following a similar strategy of interpolating and 
redshift-correcting clusters as used above. Appropriate 
cosmological dimming and redshift are then applied to 
determine the contribution of the slice to the sky ob- 
served at z = 0. We generate a radio image for each 
cluster by projecting its density and turbulent pressure 
onto the skymap and computing the relevant radio inten- 
sity using a given set of radio model parameters, ensur- 
ing that the integrated radio power across the projected 
cluster is equal to the value obtained using My and F^ 
in the above sections. We only project gas values within 
Ry. For halos not within the high-resolution sample, we 
identify the nearest high-resolution cluster in mass and 
copy that high-resolution image to the location of the 
low-resolution halo. While this procedure is admittedly 
somewhat crude, it does allow us to explore some of the 
observational consequences of these models and demon- 
strates a method of generating radio maps in the future 
using more sophisticated and realistic simulated data. 

Figure [12] shows the entire radio sky containing our 
simulated clusters at 120 arcsecond resolution assuming 
no background (i.e., a threshold sensitivity of mJy). 
This resolution best approximates the LOFAR beam at 
an average frequency of ^ 120 MHz and a longest base- 
line of L ~ 2 km. For this example we have chosen model 
parameter set Al. This map particularly highlights the 
paucity of radio halos in the universe, even at low sen- 
sitivity thresholds, but it is useful for providing a mock 
all-sky map for linking simulations to observations. 

Figure [13] highlights a region of the sky 6 degrees on 
a side at a resolution of 10 arcseconds, representing the 
high-resolution capability between 20 and 240 MHz at 
the longest baseline configuration of LOFAR. We also 
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Fig. 10. — Radio halo total counts at 1.4 GHz (solid lines) and 150 MHz (dotted lines) versus flux limit in mJy. Error bars indicate la 
Poisson uncertainties. We have identified each model with its designation from Table [2] and the portions of the model that change in the 
given model group. (S) is given in units of fiG. 
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Fig. 11. — Radio halo counts for z < 0.2 at 1.4 GHz (solid lines) and 150 MHz (dotted lines) versus flux limit in mJy. Error bars indicate 
la Poisson uncertainties. We have identifled each model with its designation from Table |2] and the portions of the model that change in 
the given model group. (B) is given in units of /iG. 



draw contour levels at varying sensitivities: 1, 10, and 
30 mJy. These sensitivities represent different configura- 
tions of the LOFAR array. At high resolution and peak 
sensitivity, we are able to clearly distinguish several sub- 
structures and features within the two radio halos, indi- 
cating that LOFAR may be able to cleanly distinguish 
various radio power models based on their dependence 
on local gas density or local turbulent pressure, which 
can have different characteristic structures in the clus- 
ter atmosphere (Figure [s]). At lower sensitivities, we 
can still distinguish features in the cluster cores, and 
early LOFAR images of nearby and bright radio halos 
may also provide useful distinguishing results. We will 
present a detailed radio morphological study, which re- 



quires knowledge of the spatial dependence of the mag- 
netic field, in a future paper. 

Figure [T4| shows the same region of the sky as above 
with a much lower resolution of of 240 arcseconds. The 
contours are the same as above. While we lose signifi- 
cant information about distant and small clusters, some 
larger clusters, such as the one shown, still show signifi- 
cant structure even at lower resolutions. We see that we 
can still identify substructure within the large cluster, 
and the effects of higher sensitivity thresholds are lim- 
ited to distant clusters and the outer regions of nearby 
objects. These results are encouraging, since they in- 
dicate that LOFAR may be able to give detailed radio 
maps of many radio halos. 
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Fig. 12. — Example radio halo all-sky map. This map assumes mjy sensitivity threshold and a resolution of 120 arcseconds. 




Fig. 13. — Example radio halo partial-sky map at 10 arcsecond 
resolution. The color scale is the logarithm of radio power in mJy. 
Contours are drawn at 1, 10, and 30 mJy levels. 

8. CONCLUSIONS 

We have introduced the first set of radio halo statistics 
derived entirely from large-scale cosmological simulation. 
Our radio power model is sufficiently broad to encompass 
many viable and more realistic models of cosmic ray gen- 
eration and synchrotron emission in clusters of galaxies. 
Our approach demonstrates the viability of using large- 
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Fig. 14. — Example radio halo partial-sky maps at 240 arcsecond 
resolution. The color scale is the logarithm of radio power in mJy. 
Contours are drawn at 1, 10, and 30 mJy levels 

scale simulation to bridge simulations and observations, 
both by deriving radio halo statistics from the simulated 
data to constrain possible radio power models and by 
producing mock radio sky maps that can be directly com- 
pared to observations. 
From our analysis we have determined that the slope 
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and normalization of the -Pi.4GHz — and Pi. 4 — re- 
lations are potentially key probes of the various models 
of cosmic ray generation. They also allow us to place 
limits on the average cluster magnetic field strength and 
the scaling of magnetic fields with cluster mass. With 
the uncertainties from only 131 observed radio halos we 
can significantly constrain the scaling of radio power with 
cluster mass and turbulent pressure. With the 131 ob- 
jects of our high-resolution sample we are able to clearly 
separate some models with strong statistical significance. 
The evolution with redshift of these relations also allows 
us to potentially distinguish various models. Future low- 
frequency missions, which will surveys large portions of 
the sky, may then potentially capture enough objects to 
perform a similar analysis and obtain these constraints. 

Our limited volume prevents us from counting the very 
largest — and hence rarest — clusters. However, we 
predict an order of magnitude fewer high-frequency ra- 
dio halos at low mass (where the undercounting effects 
of our simulation volu me are negl i gible) than the anal- 
ysis of CSB06 and Cassano et al. (2010). Some of this 
discrepancy might be due to our lack of steep-spectrum 
halos, which get counted via inclusion of the calculation 
of the synchrotron break frequency, Vb- Instead, we just 
assign a radio halo to 5% of our clusters. However, we 
found that adjusting the spect ral index to 1.9 (i.e ., the 



average spectral index found by Cassano et al.|2010 1 only 
increased the number counts by roughly 507o, which is 
not nearly enough to explain the differences. The analyt- 
ical estimates of CSB06 may be over-estimating the level 
of turbulence in the intracluster medium, whereas our 
simulations provide direct access to the cluster dynam- 
ical state. On the other hand, our approach of refining 
only cells but not dark matter particles may somewhat 
limit our dynamic range in the clusters, leading to an 
underestimate of the turbulence (although we agree with 
other, better-resolved simulations). This issue will be 
more carefully examined in subsequent work. 

We have uncovered many degeneracies among the scal- 
ings of radio power with cluster mass and turbulent pres- 
sure and the mass-dependence of cluster magnetic fields. 
These degeneracies can be broken by several methods. 
For example, a better understanding of the relationship 
between cluster mass and magnetic field will constrain 
our (B) and b parameters, allowing us to make more con- 
clusive statements about the observational limits placed 
upon the a and c parameters. On the other hand, more 
high-resolution radio and X-ray images of clusters may 
constrain the effectiveness of the various mechanisms of 
generating cosmic rays, which would further constrain 
our a and c parameters. In particular, measurements 
of intracluster gas velocities (such as may be possible 
with future X-ray spectroscopic missions) would enable 
mapping of the projected turbulent pressure, which could 
then be compared with the projected mass distribution. 

We find that low-frequency surveys are potentially ca- 
pable of distinguishing and constraining the scaling of ra- 
dio power with cluster mass and turbulent pressure, since 
the radio halo number counts are much higher at lower 
frequencies. Even though low-frequency surveys have rel- 
atively low sensitivity, their ability to map large portions 
of the sky to moderately high redshift means that they 
can gather many more objects than high-frequency (and 
more sensitive) observations. However, our estimates in- 



dicate that LOFAR will only see on the order of 10 radio 
halos within our studied mass range, and that future mis- 
sions with more sensitivity will be required in order to 
cleanly distinguish models. Future radio observations are 
especially important since the gamma ray emission asso- 
ciated with the production of cosmic rays from hadronic 
secondary interactions from clusters might be too small 
for Fermi to detect, which means that we may n ot be able 
to us e this instrument to distinguish models (Brunetti 
20091). 



Similarly, since hadronic secondary models of cosmic 
ray production are highly degenerate with re-acceleration 
models when only considering total counts and inte- 
grated cluster quantities, high-resolution low- frequency 
radio images are required in order to effectively distin- 
guish these models. Although high-frequency observa- 
tions can also (and do) produce similar maps, surveys 
such as LOFAR have the unique capability of captur- 
ing many such images, potentially providing a statisti- 
cally significant means of distinguishing models based 
on morphological differences. Our simulated radio clus- 
ter and sky maps are freely available upon request to the 
authors or via the project website ^ under the section 
jProjects/RadioHaloMaps We have produced images 
at a variety of sky coverage areas, sensitivity limits, and 
resolutions for all of the models described above. These 
images are simple FITS files. The images are straightfor- 
ward to produce, allowing us to explore further refine- 
ments to the models and more sophisticated instrument 
modeling. 

Since this initial work is highly preliminary, we have 
room for many improvements and modifications to make 
stronger connections with observations. As an immedi- 
ate improvement we may perform simulations with larger 
volumes than our 1 Gpc box (or, equivalently, per- 
form multiple realizations of the same volume) in order 
to capture more massive objects. With larger volumes 
we can also capture more low- and moderate-mass ob- 
jects to obtain better statistics for the Pi.4ghz — My and 
other relations. These simulations will allow us employ 
more sophisticated models of cosmic ray generation and 
evolution as well as enable us to calculate statistics of ra- 
dio halo morphologies. With more simulations, we may 
begin to investigate the dependence of radio halo counts 
on cosmological parameters. We may also begin to self- 
consistently include magnetic fields, although our results 
in this approach would be tied to a specific model of mag- 
netic field injection and growth. Similarly, we can begin 
to investigate generating and propagating cosmic rays in 
the simulation, although it is difficult to scale current 
methods to large volumes and high resolutions. Since 
the expected low-luminosity radio halo number counts 
are so sensitive to the level of turbulence in cluster at- 
mospheres, we must incorporate more careful techniques 
for estimating this. However, our results demonstrate 
an important first step in bridging simulations and ob- 
servations to more fully understand the large-scale radio 
universe. 
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